Dynamic structure factor of a Bose Einstein condensate in a ID optical lattice 
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We study the effect of a one dimensional periodic poten- 
tial on the dynamic structure factor of an interacting Bose 
Einstein condensate at zero temperature. We show that, due 
to phononic correlations, the excitation strength towards the 
first band develops a typical oscillating behaviour as a func- 
tion of the momentum transfer, and vanishes at even multiples 
of the Bragg momentum. The effects of interactions on the 
static structure factor are found to be significantly amplified 
by the presence of the optical potential. Our predictions can 
be tested in stimulated photon scattering experiments. 

When a Bose-Einstein condensate is loaded into an op- 
tical lattice, its properties change in a very strong way 
[1]. For deep potential wells, it even happens that the 
coherence of the sample is lost and one observes the tran- 
sition to the Mott-insulator phase [2,3]. Interesting phe- 
nomena occur also for low optical potential depth, for 
instance Bloch oscillations [4] and tunneling effects [5-7] 
can be investigated in this regime. In ID optical lattices 
the transition to the insulator phase is expected to take 
place for very large intensities of the optical lattice, so 
that there is a very extended range of parameters where 
the gas can be described as a fully coherent system. 

In this Letter, we study the elementary excitations of 
an interacting Bose gas in the presence of a periodic po- 
tential and discuss how these states can be excited via in- 
elastic processes using, for example, Bragg spectroscopy 
[8,9]. To this purpose we develop the formalism of the 
dynamic structure factor, a quantity directly related to 
the linear response of the system. 

We will restrict ourselves to the case of a system in the 
presence of a 1-dimensional optical potential 
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created by two counterpropagating laser beams. In 
Eq.(l) d is the lattice spacing and s is a dimensionless pa- 
rameter which denotes the intensity of the laser in units 
of the recoil energy En = q^/2m. Here qs = Tin /d is 
the Bragg momentum denoting the boundary of the first 
Brillouin zone and m is the atomic mass. The inclusion 
of an additional harmonic potential produced, for exam- 
ple, by magnetic trapping does not modify the excitation 
spectrum in a profound way, unless the wave length of 
the excitation is comparable with the size of the sample. 
Along the tranverse directions we assume uniform con- 
finement, so that the 3D Gross-Pitaevkii (GP) equation 



for the ground state order parameter ip{z) takes the ID 
form 
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where n is the average 3D density and the order parame- 
ter If is normalized according to X^^^yj = 1. As 
usual g = Anh'^a/m is the interaction coupling constant 
fixed by the scattering length a. 

The elementary excitations correspond to the solutions 
of the linearized time-dependent GP-equation and are 
described by the Bogoliubov equations 
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where tp is the ground state solution of Eq.(2) and the am- 
plitudes Ujq and Vjq satisfy the normalization condition 

I-d/2 [\'^jqi^)\'^ ~ \'^j<]i^)\'^] '^^ = 1- The solutions Ujq{z) 
and Vjq{z) are Bloch waves: Ujq{z) — ejip{iqz/h)ujq{z) 
where Ujq is periodic in space with period d and analo- 
gously for Vjq. For each value of the quasi- momentum 
q, Eqs.(3,4) provide an infinite set of solutions ^j{q), 
forming a band structure labelled with j ("Bogoliubov 
bands"). Due to the periodicity of the problem, the so- 
lutions of Eqs.(3,4) with q restricted to the first Brillouin 
zone and j varying over all the bands, exhaust the ele- 
mentary excitations of the system. Still it is often con- 
venient to consider values of q outside the first Brillouin 
zone and to treat the energy spectrum and the functions 
M,„ and periodic in quasi-momcntum space with 

period 2qB (see Fig.l). The intensity s of the optical po- 
tential and ratio gn/ Er between the interaction and the 
recoil energy are the relevant dimensionless parameters 
in terms of which we will discuss the physical behaviour 
of the system. First numerical solutions of Eqs.(3,4) in 
the presence of a periodic potential were obtained in [10]. 

The capability of the system to respond to an excita- 
tion probe transferring momentum p and energy Tilo is 



1 



described by the dynamic structure factor. In the pres- 
ence of a periodic potential the dynamic structure factor 
takes the form 



S {P, t^) = 51 (P^l^ ^'^ ~ '^3 ) ' 



(5) 



where Zj{p) are the excitation strengths relative to the 
j*^ band (sec Eq.(ll) below) and hujj{p) arc the corre- 
sponding excitation energies, defined by the solutions of 
Eqs.(3,4). Note that p, here assumed to be along the 
optical lattice {z axis), is not restricted to the first Bril- 
louin zone, being the momentum transferred by the ex- 
ternal probe. In this respect, it is important to point 
out that, while the excitation energies h.ujj{p) are peri- 
odic as a function of p, this is not true for the excitation 
strengths Zj. 

The dynamic structure factor satisfies important sTim- 
rules. The integral of the dynamic structure factor pro- 
vides the static structure factor (non energy weighted 
sum-rule) 
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As we will see later S{p) is strongly affected by the com- 
bined presence of two-body interactions and of the optical 
lattice. 

A second important sum-rule obeyed by the dynamic 
structure factor is the model independent /-sum rule 



/ 



hujS{p, u))duj 



P 

2m' 



(7) 



Another important sum-rule is the compressibility 
sum-rule (inverse energy weighted sum-rule) 
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where k = [n(dii/dn)]~^ is the thermodynamic com- 
pressibility. The density dependence of the chemical 
potential, and hence k, can be obtained by solving the 
GP equation (2) [11]. The compressibility k is naturally 
expressed in terms of the sound velocity c, characteriz- 
ing the low q phononic behaviour of the dispersion law 
{hw = cq), through the relation 



1 



K = 



(9) 



where the effective mass m* differs from the bare mass 
because the Hamiltonian is not translationally invariant. 
The effective mass in the presence of the external poten- 
tial (1) has been recently calculated in [12]. 

In a uniform Bose gas, the sum (5) is ex- 
hausted by a single mode with energy hujBip) = 
y^p2/2m (p2/2m -|- 2gn). In this case the static structure 
factor obeys the Feynman relation 



Sb{p) = 
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where we have used the /-sum rule (7). For p ^ the 
static structure factor (10) behaves like |p|/2mc_B, while 
the compressibility sum-rule (8) becomes l/2mc^, where 
cb = \J gn/m is the Bogoliubov sound velocity. The 
suppression of Sb (p) at small momenta is a direct conse- 
quence of the phononic correlations. For large momenta, 
instead, the static structure factor (10) approaches unity 
(see dotted lines in Figs. 2). Notice that in the absence of 
two-body interactions, S{p) = 1 for any value of p. 
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FIG. 1. Bogoliubov bands for s = 10 and gn = 0.5Er; 
excitation strengths Zj towards the states in the first three 
bands for p = —1.2qB, p = O-Sqa and p = 2.8gs- 

In the presence of the optical lattice the behaviour of 
the dynamic structure factor changes in a drastic way. 
In particular, for a given value of momentum transfer 
p, it is possible to excite several states, corresponding 
to different bands (see Eq.(5) and Fig.l). An important 
consequence is that on one hand it is possible to excite 
high energy states with small values of p, and on the 
other hand one can excite low energy states, belonging 
to the lowest band, also with high momenta p outside 
the first Brillouin zone. This behaviour introduces new 
possibilities in Bragg spectroscopy experiments. 

In general, the dynamic structure factor has to be 
calculated numerically. Starting from the solution of 
Eqs.(3,4), the excitation strengths Zj can be evaluated 
using the standard prescription (see for example [13]) 



d/2 



d/2 



(11) 



where q belongs to the first Brillouin zone and is fixed by 
the relation q = p ^ 2£qB with £ integer. This equation 
shows that, by solving Eqs. (3,4) within the first Bril- 
louin zone, one can calculate the strength Zj{p) also for 
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values oip outside the first Brillouin zone. In Figs.2(a,b), 
we show our results for two different choices of the inter- 
action at s = 10. 




FIG. 2. Static structure factor (full line), Zi{p) (dashed 
line) for s = 10 and static structure factor in the uniform gas 
(s = 0, dotted line); for gn = 0.02£ij (a) and gn = O.^Er 
(b). 



Let us first discuss the dynamic structure factor in the 
low energy region and in particular the behaviour of the 
contribution Zi{p) arising from the first band (dashed 
line). We find that Zi{p) exhibits characteristic oscil- 
lations whose amplitude is suppressed at large p. The 
zeros of Zi{p) at p = 2iqB {(- integer) reflect directly 
the phonon behaviour of the excitation spectrum which 
vanishes at the same values (see Fig.l). 

The behaviour of Zi{p) can be studied analytically in 
the large .s limit where the tight binding approximation 
applies. In this limit one can approximate the solutions 
of Eqs. (3,4) in the lowest band by 



««(^) = U,iY. (^""""""fi^ - kd) (12) 

k 

and analogously for Vq{z), where f{z) is a function local- 
ized near the bottom of the optical potential V a,t z = 0, 

and k labels the potential wells. Within this approxima- 
tion the function / also characterizes the ground state 
order parameter which reads (p{z) = J2k ./(-^ ^ ^d). 

In the tight binding approximation the dispersion law 
of the lowest band takes the Bogoliiibov-like form [14] 

nw{p) = Ve(p)Kp) + 2«-i], (13) 

where 

e(p) = 2 S sin^ (g) . (14) 

In the above equations S is the tunneling rate of parti- 
cles between two consecutive wells. The tunneling rate 
is related to the effective mass entering the compressibil- 
ity sum-rule (8,9) hy 6 = 2mER/'K'^m* and decreases 
by increasing the laser intensity s. The parameter k 
is the compressibility of the gas as emerges from the 
low momentum behaviour of the dispersion law (13): 
fiu) = ^ KT^ jm*p. In the tight binding limit, one finds 

K-^=9ndl%nz)dz [15]. 

By approximating the function f{z) with the gaussian 

/(z) = exp[— z^/2(T^]/(7r^/^Y^), one finds, after some 
straightforward algebra, the result 

for the strength relative to the first band, where the 
width a can be calculated numerically by minimization of 
the ground state energy and behaves like a ~ s~^l^dl'K 
for s \. Equation (15) reproduces with good accu- 
racy the numerical results obtained by solving the Bo- 
goliubov equations for relatively large values of s. It 
accounts for both the suppression of the strength at 
large p through the gaussian term, and the oscillating 
behaviour through the Feynman-like term e{p) /?iu}{p). 
The strength Zi has a maximum close to the edge of 
the first Brillouin zone, where it takes approximately the 
value Zi{qB) « \/k5/{k5 + 1). This simple expression 
shows that Zi is quenched both by increasing interac- 
tions [k 0) and by increasing the optical potential 
((5^0). In the cases of Figs.l(a,b) one has k,6 = 0.95 and 
0.056 respectively. In the noninteracting case (k~^ = 0) 
one has e{p) = huj{p) and the strength (15) reduces to 
^i{p) = exp (— 7r^cr^p^/2d^gg). The comparison then 
clearly shows that the oscillating behaviour of Zi{p) as 
well as its quenching at large s are a direct consequence of 
two-body interactions. On the other hand, two-body in- 
teractions scarcely affect the strengths towards the higher 
bands, provided gn <C ^/sEn. 
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The quantity Z\ (p) could be measured in Bragg spec- 
troscopy experiments by tuning the momentum and the 
energy transferred by the scattering photon to the values 
of J) and fiLO corresponding to the first Bogoliubov band. 
In order to detect a sizeable signal at large p and to point 
out the corresponding oscillating behaviour of Z\ (p) (see 
Figs. 2), the intensity s of the optical potential should be 
neither too small, nor too large. In fact for s ^ the 
strength to the lowest band becomes weaker and weaker 
if p is outside the first Brillouin zone. For large s the 
strength is instead quenched for all values of p because 
of the presence of two-body interactions. 

In Figs. 2 we also report the results for the static struc- 
ture factor (full line) corresponding to the sum S{p) = 
Zj{p). One finds that for weak interactions (Fig. 2(a)) 
also the static structure factor exhibits characteristic os- 
cillations, reflecting the contribution from the first band. 
This effect is less pronounced for larger values of gn 
(Fig. 2(b)) due to the quenching of Zi{p). In both cases 
one observes an important difference with respect to the 
behaviour of S{p) in the uniform gas (10) (dotted lines 
in Figs. 2). 

The behaviour of S{p) at small momenta can be 
described exactly using sum-rule arguments. In fact, 
phonons exhaust both the non energy and inverse en- 
ergy weighted sum rules when p — > 0, high energy bands 
giving rise to contributions of order p^ . As a consequence 
high energy bands contribute to the /-sum rule (7) but 
can not affect the low p behaviour of the non energy 
weighted moment (6) which behaves like |p|, nor the in- 
verse energy weighted moment (8) which approaches a 
constant value when p ^ 0. The result is that the low p 
behaviour of the static structure factor is entirely deter- 
mined by phonon correlations and behaves like 

.^0 ^ (16) 

consistently with the phononic dispersion law and 
Eqs.(8,9) for the compressibility sum-rule. It is worth 
noticing that result (16) holds for any value of s and gn. 
In the absence of the optical lattice one has m* — m and 
c coincides with the Bogoliubov sound velocity cb ■ Since 
one can write to*c = V ■m*K~^ and both m* and in- 
crease with s, one finds that the presence of the lattice 
results in a suppression of the static structure factor at 
low values of p, as clearly shown in Figs. 2. 

Let us conclude by recalling that in Bragg scatter- 
ing experiments one actually measures the imaginary 
part of the response function x rather than the dy- 
namic structure factor. The two quantities are related 
by the equation Im(x) = — 7r(S'(p, w) — S{—p, —oj)). The 
subtraction between the two terms can be crucial also 
at low temperatures. Actually, due to thermal excita- 
tion of phonons, the dynamic structure factor exhibits 
a strong temperature dependence when fiu < ksT, 
even if T is much smaller than the critical temperature 



for Bose-Einstein condensation. However the difference 
S{p, Lo) — S{—p, — w) cancels out most of this temperature 
dependence so that the measurement of Im(x) provides 
reliable information on the zero temperature behaviour 
of the dynamic structure factor [18]. 
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